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ABSTRACT 

We explore the vicinity of the Milky Way through the use of spectro-photometric data from the Sloan 
Digital Sky Survey and high-quality proper motions derived from multi-epoch positions extracted from 
the Guide Star Catalogue II database. In order to identify and characterise streams as relics of the 
Milky Way formation, we start with classifying, select, and study 2417 subdwarfs with [Fe/H] < —1.5 
up to 3 kpc away from the Sun as tracers of the local halo system. Then, through phase-space 
analysis, we find statistical evidence of five discrete kinematic overdensities among 67 of the fastest- 
moving stars, and compare them to high-resolution N-body simulations of the interaction between a 
Milky-Way like galaxy and orbiting dwarf galaxies with four representative cases of merging histories. 
The observed overdensities can be interpreted as fossil substructures consisting of streamers torn 
from their progenitors; such progenitors appear to be satellites on prograde and retrograde orbits on 
different inclinations. In particular, of the five detected overdensities, two appear to be associated, 
yelding twenty-one additional main-sequence members, with the stream of Helmi et al. (1999) that our 
analysis confirms on a high inclination prograde orbit. The three newly identified kinematic groups 
could be associated with the retrograde streams detected by Dinescu (2002) and Kepley et al. (2007); 
whatever their origin, the progenitor(s) would be on retrograde orbit(s) and inclination(s) within the 
range 10° A 60°. 

Finally, we use our simulations to investigate the impact of observational errors and compare the 
current picture to the promising prospect of highly improved data expected from the Gaia mission. 
Subject headings: Galaxy: formation — Galaxy: halo — Galaxy: kinematics and dynamics 


1. INTRODUCTION 

The formation and evolution of galaxies is one of the 
outstanding problems in astrophysics, one which can be 
profitably engaged directly through de tailed stu dy of our 
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ical structure formation, 
galaxies such as the Milky Way grow by mergers and ac¬ 
cretion of smaller systems, perhaps similar to what are 
now observed as dwarf galaxies. These satellite galaxies 
- torn apart by the tidal gravitational field of the parent 
galaxy - are progressively disrupted, giving rise to trails 
of stellar debris streams along their orbits, spatial signa¬ 
tures that eventually disappear due to dynamical mixing. 
After the accretion era ends, a spheroidal halo-like com¬ 


ponent is left from their collective assemb l y (e.g., Searle 
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Of all the Galactic components, it is indeed the stellar 
halo that offers the best opportunity for probing d etails 
of the merging history of the Milky Way (see, e.g., Helmi 


2008). Past explorations have demonstrated that there 


is a concrete possibility to identify groups of hal o stars 
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much like those seen in the halo system of the Milky Way, 
are sensitive to recent (within the last 8 Gyr) merging 
events, and are more prominent in the outer region of the 
halo (galactocentric radii beyond 15 — 20 kpc), whereas 
the inner-halo region appears significantly smoother. 
Based on data from the SEGUE spectroscopic sur- 
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main-sequence turnoff stars in the inner-halo region of 
the Milky Way (within ~ 20 kpc from the Sun) exhibit 
clear evidence for radial velocity clustering on small spa¬ 
tial scales (they refer to these as ECHOS, for Elements 
of Cold HalO Substructure). They estimated that about 
10% of the inner-halo turnoff stars belong to ECHOS, 
and inferred the existenc e of about 1000 ECHOS i n the 
entire inner halo volume. |Schlaufman et al. (2011) sug¬ 


gest that the most likely progenitors of ECHOS are dwarf 
spheroidal galaxies with masses on the order of 10 9 M 0 . 

In the Solar Neighbourhood, up to 1 — 2 kpc of the Sun, 
stellar streams have also been discovered as overdensities 
in the phase-space distribution of st ars, i ntegrals o f mo¬ 


tion and action-angle variables (see |Klement 12010", for 


a review). Prominent examples are the two stellar de¬ 
bris streams in t he halo population pa ssing close to the 


Sun detected by Helmi et al. (1999) when combining 


high-quality HIPPARCOS proper motions with ground- 
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based observations. Formed via destruction of a satellite 
whose debris now occupy the inner halo region with no 
apparent spatial structure, these streamers retain very 
similar velocities and are seen as clumps in angular mo¬ 
mentum space where s tars from a common progenitor 
appear rather confined (Helmi fe de Zeeuw] 2000). 
Besides the Helmi stream, u> Centauri (Bincscu ]|2002 


Ma.jewski et al. 2012), the Kapteyn an d Arcturus 


(e.g., Eggen 197ip streams, Klement (2010) lists a few 


other halo substructures found in the solar neighborhood: 
these are still small numbers compare d to the few hum 
dred streams expec ted (i.e. 300-500, |Helmi & White 


1999 Gould 


only be broken with 6D (phase-space) or 7D (includ¬ 
ing abundances) information achievable by integrating 
astrometry, photometry, and spectroscopy. 

The SDSS - GSC II Kinematic Survey (from now on 
SGK S) we exploit here wa s produced to serve this task 


(see Spagna et al. 2010bI. In the future, ne w ground- 


and space-based surveys such as G aia (e.g., Perryman 


et al. nj|001|_^ronet^ 


GFS; Gilmore et al 


. || 2012 |, a 

leh-precisi 


2005), Gaia ESU Survey 
and UAMOST (Zhao et al. 


2012) will provide high-precision data that will usher us 


in a new era of Milky Way studies. 

In Sect. [2j we introduce the data used to isolate a sam¬ 
ple of nearby halo subdwarfs from the SGKS catalogue. 
The kinematic and orbital properties of the local halo 
subdwarf population are discussed in Sect. [3j where we 
present algorithms to searc h for kinematic substructures, 


over densities. In Sect. 


recovering known streams (Helmi et al. 
2002] |Kepley et al. 2007^ 
sect. [1J we ] 


1999 


Dinescu 


as well as new kinematic 


present the high-resolution 
N-body numerical simulations of four minor mergers used 
to study galaxy interactions and the properties of accre¬ 
tion events in the vicinity of the Sun. In Sect. [5] we 
investigate the impact of observational errors resulting 
from current ground-based data and from high accurate 
data expected from the Gaia satellite. Finally, in Sect. [6j 
we compare observations to these numerical simulations 
and infer the nature of the detected fast moving groups. 

2. THE SDSS - GSC II KINEMATIC SURVEY 
(SGKS) 

This study is based on a new kinematic catalogue, 
derived by assembling spectro-photometric stellar data 
from the Seventh Dat a Release of the Sloan Digital Sky 


Survey (SDSS DR7; Abazajian et al. 2009), which 


included data from the Sloan Extension ior Galactic 


Un ders tanding and Exploration (SEGUE; Yanny et al. 

n [ 2009 ] ) 


supplemented by astrometric parameters ex¬ 
tracted from the database used for the construction of 


the S econd Guide Star Catalogue (GSC II; Lasker et al. 


2008). This SDSS - GSC II catalogue contains positions, 


proper motions, classification, and ugriz photometry for 
77 million sources down tor~ 20, over 9000 square de¬ 
grees. 

Proper motions are computed by combining multi¬ 
epoch positions from SDSS DR7 and the GSC II 
database; typically, 5 — 10 observations are available 
for each source, spanning up to 50 years. The typical 
formal errors on those proper motions are in the range 
2 — 3 mas yr -1 per coordinate for 16 < r < 18.5, compa¬ 


rable with the intern al precisio n of the SPSS proper mo¬ 
tions computed by Munn et al. (2008). Although much 
of the photographic material (Schmidt plates) u sed to 
derive fi rst ep och information is in common with |Munn| 


2003|). Actually, recovering fossil struc- et al. 
tures in the inner halo is considerably more difficult, as 
strong phase-mixing takes place. This degeneracy can 


et al. (2008), plate digitisation and measurement pro¬ 
cesses, alnd the calibration methods that led to the first 
epoch positions were somewhat different, with the Munn 
et al.’ data coming from the USNO-B project. Of par¬ 
ticular relevance is the minimisation of systematic errors 
that can affect proper motion accuracy, a true driver in 
analysis like those conducted in this study. An accu rate 
validat ion of our proper motions was discussed in|Spagna| 


(2010a). 

Radial velocities and astrophysical parameters are 
available for about 151 000 sources cross-matched with 


the SDSS spectroscopic catalogue Typical accuracy are 
of 5 — 10 km s -1 in line of sight velocity, 250 K in effec¬ 
tive temperature, T e g, 0.25 dex in surface gravity, log g , 
and 0.20 dex in metallicity, [Fe/H], as estimated within 
the SEGUE Spectral Parameter Pipeline (SSPP; i. e.. 
Re Fiorentin et al. 2007| |Lee et al. 1 2008a|b Allende 


Prieto et al. |2008p . We specify that the sample includes 


only objects with no problems related to the spectrum 
and classified without any cautionary flag by the SSPP. 
In case of multiple spectra, we take the spectrum with 
the highest signal-to-noise ratio. 

From the SDSS - GSC II catalogue, we select as tracers 
sources with 4500 K < T e g < 7500 K and log g > 3.5, 
corresponding to FGK dwarfs. 

The observed magnitudes are corrected fo r interste llar 
absorp tion via the extinction maps ofJSchlegcl e t al. 
(1998) based on the 6.1' resolution COBE/DIRBE dust 
map, that we preferred to the more recent, but of inferior 


resolution (7'- 14'), reddening maps published by Schlafly 
et al. (2014). Then, we transformed the E(B — V) to 


the SDSS photometric system by adoptin g the extinction 


ratio A r /Ay = 0.875 (from Table 1 of Girardi et al. 


20041, that is appropriate for our FGK dwarf sample. 

Photometric distances good to crd/d ~ 20% are com¬ 
puted by means of the photometric parallax rela tion es- 
tablished for FGK main-sequence stars by I vezi c et a l.] 
(2008). Here, the metallicity-dependent absolute mag¬ 
nitude relations, M r = f(g — i, [Fe/H]), use the spec¬ 
troscopic [Fe/H] instead of the photometric metallicity 
adopted by Ivezic et al. ( 2008|) . We also apply the ad- 


ditional colour thresholds from (Klement et al. (2009) in 
order to remove turn-off stars, whose estimated M r may 
be affected by residual systematic errors. 

Galactic space-velocity component^] are estimated un¬ 
der the assumption that the Sun is at a distance of 8 kpc 
from the centre of the Milky Way, the Local Standard 
of Rest (LSR) rotates at 220 km s” 1 about the Galactic 
center, and the peculiar velocity of the Sun relative to the 
LSR is (U V, Wy = (10.00,5.25,7.17) km s" 1 (Dehnen 


& Binney 1998). 

Finally, in order to minimise the effect of outliers (e.g. 
mismatches, blends and sources with low S/N) and there¬ 
fore obtain a sample with accurate distance and kine¬ 
matics suitable for our stellar stream search, we impose 

1 Throughout this work, U, V, and W indicate Galactic velocity 
components relative to the Local Standard of Rest and follow the 
convention with U positive toward the Galactic center, V positive 
in the direction of Galactic rotation, and W > 0 towards the North 
Galactic Pole. 
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TABLE 1 

Halo velocity parameters. 


(U) 

(V + 220) 

(W) 

&U 

(T V 

aw 

PUV 

Puw 

Pvw 

(km s - 1 ) 

(km s -1 ) 

(km s 1 ) 

(km s -1 ) 

(km s -1 ) 

(km s -1 ) 




15 + 2 

25 + 2 

-4 + 2 

126+ 1 

100+ 1 

91 ± 1 

-0.09 + 0.02 

-0.18 + 0.02 

0.05 ± 0.02 


Note. — The Milky Way halo velocity parameters as determined from our selected sample of 2417 FGK subd¬ 
warfs. The table lists mean velocities, dispersions, and corresponding correlation coefficients in Galactic coordinates. 
Noticeable is the correlation between U and W (see text). 


a threshold on proper-motion errors (< 10 mas yr -1 per 
component), constrain magnitudes to the range 13.5 < 
g < 20.5, limit the errors on the derived velocity compo¬ 
nents to better than 50 km s -1 , and remove total space 
velocities above 600 km s -1 . These are the properties of 
the 24 634 stars listed in the SGKS catalogue. 

3. DATA ANALYSIS 

Among the full sample of FGK dwarfs from the SGKS 
catalogue, we have selected specific sub-samples of trac¬ 
ers of the Galactic halo population in the inner-halo re¬ 
gion, and analysed their phase-space distribution. 

Here, we focus on a sample of 2417 metal-poor stars 
([Fe/H] < —1.5) outside the Galactic plane (|z| > 1 kpc), 
and located within 3 kpc of the Sun. Within this volume, 
the selected sample has median errors on U, V, and 
W of 12, 13, and 9 km s -1 , respectively; this results 
in errors in the velocity difference between stellar pairs 
not exceeding ~ 20 km s -1 . Such a value is suited for 
careful investigations of substructure, as the kinematic 
analysis presented below will show. 


3.1. Local halo velocity distribution 

From the selected sample we measure the mean ve¬ 
locities (( U),(V + 220), (W)), the velocity ellipsoid 
(cry, ay, and the correlations among velocity com¬ 
ponents (puv , Puw j Pvw) as reported in Table[l] 

The kinematic properties of the selected tracers are 
representative of the halo populati on in the vicinity of 


the Sun (e.g., Chiba & Beers 


2000 ). 


The significant correlation puw = —0.18 ± 0.02 be¬ 
tween the radial and vertical velocity components indi¬ 
cates a tilt of the velocity ellipsoid (Fig. [ T| ri ght panel). 

Usi ng the tilt formula (see, e.g., Binney & Merrifield 
19981 


tan2<5j/w = 


uw 


2 puw &u&w 
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and the values in Table [l] for the correlation coefficient 
and velocity dispersions along the U and W axes, we 
derive a tilt angle of Suw = —14.5° ± 1.4°, revealing 
that the ((7, W) distribution points toward the Galactic 
center. 

In fact, for our halo sample of 2417 FGK subdwarfs 
with (z) ~ 1.2 kpc and (R) ~ 8.3 kpc, we estimate a 
mean position angle (tan ~ 1 (z/R)} ~ 8.3°. This result 
is fairly consistent with the tilting effects on the velocity 
ellipsoids due to the gravitational potential prod uced by 
the stellar disk and dark matter halo (Bond et al. |2010, 
and references therein). We also measure smaller but 


stastistically significant correlations in (U, V) and (U, W) 
velocity-planes. 

In the following, we look for halo streamers in the high 
velocity tail of the ( U , V, W) velocity distribution, where 
kinematic substructures are more easily detected. In or¬ 
der to select high velocity stars, we model the velocity 
distribution as a tilted Schwarzschild ellipsoid: 

f(U,V,W) = const • e ~i E ^ u,v,w ^ (2) 

where E is the velocity function defined by: 

E(U,V,W) = 


Ruu 
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Z R 


( U-(U) \ 2 I «vv ( V-{V) \ 2 | R ww ( W-{W) \ 2 
\ a v ) ^ R \ a v ) ^ R \ a w ) 

■( 


U-(U) 

cru 


V-(V) 

a v 


/ E-(V) \ ( W-(W) \ 

V °v J V aw ) 


o Ruw 
Z R 


U-(U) 
. a v 


W-(W) 

a w 


( 3 ) 

Here, R represents the determinant of the symmetri¬ 
cal matrix 1Z of the correlation coefficients pij = Rij /R 
(for i,j = U, V + 220, W), and Rij designate the cofac- 
tor of t h e corresp on ding c orrelation element in 7 Z (e.g., 
Trump ler fc Weaver 1 1953]). 

Figure |T] shows the kinematic distribution for the indi¬ 
vidual components, (( 7 , V +220, W), of the space-velocity 
vector for the full sample of 2417 selected halo stars; the 
242 objects comprising the sample of the 10% highest 
velocity tail are represented with crosses. 

As expected, the overall velocity distribution is rela¬ 
tively smooth, because of the strong phase-mixing that 
takes place i n the inner-halo region, and slowly prograde 
(e.g., Helmi |2008). 

However, as their motions (in direction and speed) are 
well separated from those of the other nearby subdwarfs, 
we intend to study the degree of dumpiness of the 10% 
fastest-moving objects. The case study is that, of all the 
objects passing within a few kiloparsecs of the Sun, some 
are part of a diffuse local stellar halo, while some could be 
debris of accretion events and remnants from the outer- 
halo population currently in the Solar Neighbourhood. 

Before starting to look for kinematic substructures, we 
check for thick disk stars that possibly contaminate our 
halo trac ers. Here, we a pplied to kinematic method de¬ 


scribed in Spagna et al. (2004) and estimate the fraction 
of subdwarfs that is consistent with the 3D velocity dis¬ 
tribution of the thick disk population. By assuming a 
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Fig. 1.— Distribution of nearby halo stars in velocity space for our selected sample of 2417 FGK subdwarfs, with [Fe/H] < —1.5 and 
\z\ > 1 kpc within 3 kpc of the Sun. The 10% fastest-moving stars (242) are marked as crosses. 


velocity ellipsoid, as estimated by |Pasetto et al. (2012), 
an d a rotation velocity 1 4 = 150 km s -i , as measured 
by |Spagna et al. (j2010a) for metal-poor thick disk stars 
with [Fe/HJ ~ — 1 dex, we found at a 2cr (i.e. 87%) con¬ 
fidence level a ~ 10% maximum contamination of thick 
disk in the whole sample of 2417 halo tracers. Instead, 
no contaminant is expected among the subsample of the 
10% fastest objects. 

We use the samples described above to detect and sub¬ 
sequently identify kinematic halo substructures in the So¬ 
lar Neighbourhood as groups of stars moving with similar 
velocities and directions. Detection is accomplished by 
performing a statistical test based on individual kine¬ 
matics aimed at quantifying possible deviations from a 
smooth distribution of the background halo; cluster anal¬ 
ysis in velocity space is then applied for final confirmation 
of the substructures. 

3.2. The two-point correlation function: finding the 

clumps 

The amount of kinematic substructures that cosmol¬ 
ogy might leave in the volume is quantified by means of 
the cumulative two-point correlation function, £(v), on 
the paired velocity difference v = v, : — v ,■ that measures 
the excess in the number of stellar pairs moving within 
a given velocity difference when com pared to a represen¬ 
t ative random smooth sample (cfr. Re Fiorentin et ah 


2005 for more details). Here the random points were 


drawn from a multivariate distribution obtained from the 
observed data set by random permutations of the order 
of the velocity components V + 220 and W, after fixing 
U ; finally, the actual random pairs are obtained after 
averaging over ten independent realisations. 

This function is computed over the full sample of 2417 
halo stars, and separately for the sub-sample of the 242 
fastest-moving stars, corresponding to the 10% high- 
velocity tail. 

A statistical excess of stars with small pairwise veloc¬ 
ity differences indicates the presence of likely streamers 
made of objects with coherent kinematics. 

Figure[2]shows, using bins of 10 kms' 1 widtlFJthe two- 
point correlation function £(v) for the full sample of 2417 

We fixed the bin width following the rule that the interval 


halo stars (dots) and for the subset of the 10% fastest- 
moving stars (diamonds). While weak for the full sample, 
there is a statistically significant signal (SNR > 4) for 
the subset of the fastest stars, that peaks at 40 kms -1 : 
the excess of pairs of stars with similar velocities is very 
noticeable, and is a direct indication of the presence of 
kinematical clumps. 



IU; - 


(km/s) 


Fig. 2.— Cumulative velocity correlation function for the full 
sample of halo stars (dots), and the 10% fastest-moving subset (di¬ 
amonds) shown in Fig. [T] The error bars are derived from Poisson’s 
statistics of the counts. 


In the following, among the sample of the 242 fastest 
stars, we focus on the objects with paired velocity differ¬ 
ences less than 40 kms -1 , which yield the statistically 
significant signal seen in Fig. [21 In addition, we exclude 
isolated pairs, i.e., “groups’^ with only two objects. 
This further selection certainly reduces the number 
of detected members, however it makes the following 
analysis more robust by decreasing the contamination 
of false positives. The final sample is made of 67 stars. 


sampled is divided into as many bins as the square rooth of the 
sample size, in our case ~ 400 kms -1 /\/2417 ~ 10 kms -1 . 
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Fig. 3.— Distribution of the high-velocity tail from our selected sample with [Fe/H] < —1.5 and |z| > 1 kpc within 3 kpc of the Sun (see 
Fig. 0 Shown are 242 objects, the 10% fastest-moving stars. Among them, star symbols identify the 67 sources with pairwise velocity 
differences below 40 km s — , the stars belonging to isolated pairs have been excluded. Different colours are used to indicate stars associated 
with the five clumps recovered by the clustering analysis. 


3.3. Clustering analysis: assigning membership 

In order to classify these 67 objects, we perforin 
K —medoids clusteringj+n the 3D velocity space that de¬ 
fines the number of kinematic substructures and their 
members. This unsupervised learning algorithm is able 
to group data into a pre-specified number of clusters that 
minimises the RMS of the distance (in velocity space) to 
the center of each cluster. 

The original data set is initially partitioned into clus¬ 
ters around K data points referred to as the medoids, 
then an iterative scheme (PAM, for Partitioning Around 
Medoids) is applied to locate the medoids that achieve 
the lowest configuration “cost”. The algorithm employed 
by PAM, similar to the A'-means clustering algorithm, is 
more robust to outliers and obtains a unique partitioning 
of the data without the need for explicit multi ple startin g 
points for the proposed clusters (see, e.g., Kaufmann & 
Rousseeuw ||1990 Hastie et al. |2001 ). 

There is no general theoretical solution for finding the 
optimal number of clusters for any given data set. In¬ 
creasing K results in the error function values formally 
much smaller, but this increases the risk of overfitting. 
In order to keep the final identification safe and simple, 
we compared the results of runs with different K classes, 
and the best K resulted following visual inspection of 
the generated distribution. The solution adopted here 
is with I\ = 5: for K < 5 the clusters returned by the 
algorithm would contain a mixture of the natural under¬ 
lying groups (e.g., prograde and retrograde members in 
the same kinematic clump); for K > 5 natural groups 
further partition into “artificial” subgroups. 

The five kinematic substructures detected are visu¬ 
alised in Fig. [3] with different colours. 

The individual kinematic properties of the 67 stars be¬ 
longing to the five kinematics groups are listed in Table[2 

Other methods have been utilised to isolate groups o: 
stars in the halo like, e.g., the interesting approach specif- 
ically developed for the Virgo stellar stream by |Duffau 


et al. (2014|. On the other hand, the fact that we have 
the complete set of 3D kinematical data and that the 
whole sample is confined within 3 kpc from the Sun (i.e., 
the distance segregation is implicitly implemented in our 
sample) suggests the direct use of a classical clustering 
algorithm like PAM as the method of choice. 

3.4. Angular momentum and orbital properties 

The space of adiabatic invariants allows better identifi¬ 
cation of the different possible merging events that might 
have given rise to the observed substructures. Clump¬ 
ing should be even stronger, since all stars originating 
from the same progenitor should have very similar inte¬ 
grals of motion, resulting in a superposition of the corre¬ 
sponding streams; that is, the initial clumping of satel¬ 
lites are present e ven after the system has completely 
phase-mixed (e.g., |Helmi & de Zeeuw ||2000|). 

In this study we focus on the plane denned by the 
components of angular momentui+Jin and out the plane 
of the Galaxy’s disk, i.e. L xy and L z respectively. 

Since for a local sample ( x,y,z ) ~ (r Q ,0, 0), we re¬ 
mark that L xy ~ fq|u 2 | is dominated by the velocity 
perpendicular to the plane and L z ~ rQV y measures the 
amount of rotation of a given stellar orbit. Essentially, 
stars with high/low L xy are on high/low-inclination or¬ 
bits; stars with L z < 0 are on retrograde orbits and stars 
with L z > 0 are on prograde orbits. 

Figure 0] shows the distribution of the selected sample 
within 3 kpc of the Sun in the angular momentum dia¬ 
gram L z versus L xy . As in Fig.[3| the 10% fastest-moving 
objects are plotted as crosses, and with the sta r sym bols 
we mark the group members identified in Sect. [372] Dif¬ 
ferent colours indicate the stars associated with the dif¬ 
ferent lumps recovered by the cluster analysis in velocity 
space. 

The solid lines show the loci of the known kinematic 
structures detected by 


O 

J We used the implementation of the K —medoids c luste ring de¬ 

veloped as part of the H Project for Statistical Computing: www.r- 
projcct.org 


4 Remind that: L x = yv z — zv y , L y = zv x — xv z , L z = xv y — yv x , 
and L xy = L% + L^. Here v x = —U, v y = V + 220, and v z = W. 
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Fig. 4.— Distribution of the selected sample of 2417 FGK subdwarfs within 3 kpc of the Sun in the space of adiabatic invariants. 
Cross symbols indicate the 10% fastest-moving objects. Among them, star symbols identify the 67 sources with paired velocity differences 
below 40 km s -1 . As in Fig. |3j different colours are used to indicate stars associated with the five clumps recovered by the clustering 
analysis in velocity space. At L z > 0, the solid box shows the locus of the halo stream discovered by Helmi et al. |(|1999): the kinematic 
substructures, pink and blue stars (i.e., Groups 1 and 2 in Table |2|> , on prograde orbits in deed cover the same region. At L z < 0, the solid 
box at L z ~ —2375 shows the locus of the substructure detected oy Kepley et a l. | (|2007|) as retr ograde outliers, while the solid contour at 
L z rsj —400 identifies the c o Cen substructures remapped from the L z -L region m Dinescu ](|2002|). The area within the dashed line includes 
the kinematic Groups 3, 4, and 5 of Table [2] represented by the red, green, and yellow stars, respectively. 


Helmi et al. 

(1999 

) at 

300 < L z < 1500 kpc km s 1 and 
1400 < L xy < 2500 kpc km s -1 , 

Kepley et al. ( 

2007) at 

-3000 < L z < 
0 < L xy < 250C 

— 1750 kpc km s -1 
kpc km s^ 1 , 

)afl 

Dinescu 

to 

o 

o 

to 


-600 < L 
0 < L 


and 


< —200 kpc km s 
kpc km s — . 


and 


/ r Ltm 
^ -L'xy 


The most noticeable feature in Fig. [4] is certainly the 
kine matic g rou p corr esponding to the stream found by 


Helmi et al. (1999). Here, we ident ify 25 su bd warfs , 
includ ing 4 stars already detected by |Klement et _al. 


(2009), and 21 new members. By inspection o:FFig.|3]we 


notice that the 10 members belonging to Group 1 (pink 
star symbols) run along near-parallel orbits and cross the 


5 For the [Dinescu | (j2002| region, the curve delimiting the L xy 
upper part was derived by remapping the L,-L region shown in 
Fig. 4 of that paper; therefore 671 < L^IJ 1 < 922. 


Milky Way’s disk at high speed from South to North, and 
the 15 objects in Group 2 (blue star symbols) cross the 
Milky Way’s disk at similar speed and angle, but from 
North to South. 

The three remaining lumps of fast-moving stars (red 
Group 3, green Group 4, and yellow Group 5) appear 
on the retrograde side of Fig. [4] The pentagonal box 
confined by the dashed line includes most of the members 
of Group 3 (5 stars), Group 4 (21 stars), and Group 5 
(16 stars). 

These groups, and in particular the small Group 3, do 
not appear to be easily associated with known streams 
and, in Sect. [5] we discuss the possibility that all these 
stars come from a common progenitor or from three dif¬ 
ferent merging events. 

Anyhow, we note that Group 4 might be the par¬ 
ent populations of the counter- rotating “outliers” , with 


Vrf, < —250 km s , found by Kepley et al. (2007) 


while the slightly retrograde Group 5 (V^ ss —50 km s 
seems to be rela ted to the kinematic struc ture found by 
Dinescu (2002 |), and confir med by both Meza et al. 


(2005j 


and 


IVlajewski et al. (2012). These authors have 































































Inner Halo Streamers in the Era of Gaia 


7 


TABLE 2 

Main individual characteristics of the 67 fastest-moving stars found members of 5 different 

KINEMATICS GROUPS. 


Group 

ID 

U 

V + 220 

W 

L/xy 

L z 

[Fe/H] 



(km s _1 ) 

(km s _1 ) 

(km s -1 ) 

(kpc km s _1 ) 

(kpc km s -1 ) 

(dex) 

1“ / pink 

52209-0694-596 

143 ± 14 

33 ± 12 

293 ±6 

2289 

389 

-2.01 

52721-1050-418 

145 ± 7 

173 ±6 

155 ±4 

1390 

1473 

-2.02 


53315-1907-393 

170 ± 31 

54 ± 31 

285 ± 16 

2232 

614 

-1.75 


53349-2066-511 

154 ± 23 

158 ± 28 

169 ± 27 

1459 

1678 

-1.85 


54175-2472-370 

139 ± 13 

76 ± 16 

246 ±4 

2336 

634 

-1.75 


54574-2904-114 

243 ± 9 

35 ±9 

254 ± 3 

2301 

316 

-1.94 


54577-2906-088 

225 ± 12 

8± 12 

272 ±5 

2532 

198 

-2.16 


54621-2191-177 

231 ± 11 

16 ± 11 

257 ± 10 

2045 

405 

-1.81 


54624-2189-236 

169 ± 17 

84 ± 26 

261 ± 23 

1820 

590 

-1.73 


54629-2902-237 

172 ± 7 

172 ±9 

138 ±6 

1216 

1203 

-2.64 

2 a / blue 

52316-0559-336 

189 ± 3 

87 ± 11 

-310 ±5 

2545 

853 

-2.52 

53084-1368-399 

195 ± 5 

55 ±6 

-255 ± 4 

1906 

517 

-2.16 


53242-1896-109 c 

41 ± 14 

149 ± 11 

-256 ± 7 

2220 

1299 

-1.66 


53262-1900-359 

62 ± 19 

82 ± 21 

-283 ± 17 

2215 

706 

-2.09 


53293-1906-633 c 

27 ± 6 

139 ± 7 

-261 ± 5 

2281 

1210 

-1.91 


53467-2110-134 

69 ± 10 

102 ± 11 

-260 ± 3 

1921 

815 

-1.55 


53712-2314-639 

75 ±5 

103 ±6 

-265 ± 5 

2387 

936 

-2.57 


53726-2306-188 c 

55 ±9 

136 ± 10 

-283 ± 7 

2555 

1218 

-1.74 


53907-2209-540 c 

70 ± 13 

144 ± 12 

-256 ± 13 

1663 

1011 

-2.30 


54178-2452-540 

70 ± 6 

184 ± 7 

-216 ±3 

1735 

1556 

-1.51 


54380-2323-448 

47 ± 15 

170 ± 13 

-243 ± 13 

1858 

1324 

-2.12 


54479-2867-531 

33 ± 14 

129 ± 11 

-246 ± 12 

2153 

1097 

-1.64 


54530-2889-458 

3 ± 4 

148 ±9 

-247 ± 4 

2238 

1342 

-1.56 


54554-2918-615 

198 ± 15 

80 ± 15 

-276 ± 3 

2004 

604 

-1.98 


54580-2905-169 

80 ± 10 

115 ± 11 

-271 ± 4 

1923 

873 

-1.71 

3 b / red 

54179-2567-458 

-57 ± 12 

-112 ± 16 

214 ± 8 

1922 

-1076 

-1.64 

54463-2856-571 

-129 ±9 

-121 ± 22 

209 ±6 

1617 

-1089 

-2.09 


54551-2394-228 

-104 ± 10 

-130 ± 14 

230 ±8 

1924 

-1234 

-1.72 


54562-2920-596 

-81 ± 19 

-152 ± 19 

232 ± 10 

1521 

-1032 

-2.36 


54569-2900-046 

-46 ± 7 

-139 ±8 

218 ±3 

1712 

-1139 

-1.68 

4 b / green 

52059-0597-072 

125 ±4 

-157 ± 10 

-210 ±6 

1672 

-1280 

-2.09 

52338-0788-070 

180 ± 9 

-237 ± 11 

-142 ±9 

921 

-1767 

-2.07 


52942-1509-488 

179 ± 9 

-259 ± 22 

-105 ±8 

1137 

-2235 

-1.64 


53710-2310-141 

207 ± 39 

-159 ± 16 

-95 ± 20 

1195 

-877 

-1.60 


53762-2381-588 

213 ± 16 

-244 ± 29 

— 117 ± 14 

916 

-2517 

-1.85 


53800-2383-625 

172 ± 7 

-190 ± 14 

-164 ±8 

1275 

-1652 

-1.69 


53823-2240-213 

310 ±22 

-190 ± 26 

-100 ±4 

557 

-1513 

-1.70 


53874-2173-414 

279 ± 16 

-151 ± 18 

-89 ± 12 

199 

-688 

-1.57 


54154-2701-341 

171 ± 10 

-216 ± 16 

-191 ± 21 

1826 

-2236 

-1.74 


54169-2413-090 

149 ± 12 

-172 ± 12 

-226 ± 7 

1711 

-1626 

-1.77 


54234-2663-321 

90 ± 11 

-151 ± 13 

-218 ±3 

1581 

-1187 

-2.16 


54539-2894-314 

162 ±3 

-244 ± 12 

-112 ±5 

793 

-2008 

-2.12 


54539-2894-632 

305 ±6 

-168 ±6 

-99 ±3 

478 

-1299 

-1.59 


54544-2459-072 

239 ± 7 

-100 ±8 

-154 ± 7 

802 

-532 

-1.82 


54557-2177-009 

277 ± 16 

-100 ± 18 

-130 ± 10 

507 

-454 

-2.73 


54568-2899-316 

182 ± 17 

-160 ± 25 

— 117 ± 4 

665 

-1294 

-1.92 


54595-2932-091 

268 ± 17 

-120 ± 15 

-156 ±7 

737 

-704 

-1.63 


54597-2561-326 

265 ± 15 

-89 ±3 

-127 ±4 

750 

-289 

-1.96 


54616-2460-420 

186 ±8 

-169 ±9 

-169 ±9 

1180 

-1246 

-1.93 


54616-2460-616 

192 ± 16 

-199 ± 15 

-152 ± 16 

977 

-1450 

-1.98 


54631-2911-151 

182 ± 11 

-187 ± 11 

-132 ±8 

710 

-1212 

-1.96 

5 b / yellow 

53035-1433-600 

322 ± 10 

—8 ± 13 

22 ±6 

1004 

-56 

-1.58 

53240-1894-079 

277 ± 22 

-99 ± 18 

59 ± 18 

120 

-273 

-2.10 


53315-1907-353 

352 ± 35 

-33 ± 35 

68 ± 18 

258 

-32 

-1.50 


53770-2387-010 

322 ±6 

—2 ± 12 

21 ±4 

559 

-66 

-1.67 


53848-2437-060 

314 ± 13 

17 ± 22 

121 ± 13 

1787 

-228 

-2.18 


53876-2134-516 

318 ± 13 

-94 ± 13 

-17 ±5 

528 

-518 

-1.52 


53918-2539-196 

240 ± 10 

-13 ±8 

126 ±8 

1356 

121 

-1.93 


54082-2325-126 

278 ± 23 

0± 16 

115 ± 15 

468 

420 

-1.61 


54156-2393-459 

336 ± 11 

-115 ±9 

-8 ±5 

527 

-1231 

-1.74 


54243-2176-476 

350 ± 6 

-32 ±4 

50 ±4 

795 

131 

-1.73 


54271-2449-590 

261 ± 11 

-12 ±9 

147 ± 7 

1686 

321 

-2.21 


54368-2804-351 

302 ± 15 

-100 ± 19 

24 ± 19 

332 

-523 

-1.67 


54536-2871-426 

286 ± 21 

-97 ± 40 

42 ± 11 

1155 

-888 

-1.93 


54594-2965-227 

288 ± 14 

-107 ± 16 

-31 ± 6 

387 

-974 

-1.64 


54594-2965-272 

328 ± 14 

-37 ± 12 

25 ±5 

782 

-448 

-1.54 


54616-2929-342 

310 ± 7 

-58 ±8 

36 ± 2 

625 

-403 

-1.60 


Note. -— Main individual characteristics of the 67 selected fastest-moving stars with paired velocity differences 
less than 40kms —1 , and belonging to 5 different kinematics groups. T hese are 25 subdw arfs (21 new) on 2 
streamers (Group 1 and Group 2) both members of the stream known to|Helmi et al. | (|1999[) originally made of 
red giants and RR Lyrae. The remaining 42 are subdwarfs members of three newly discovered kinematic groups 
(Group 3, 4, and 5); see text for their dynamical interpretation. 

a Subdwarf members associated with the |Helmi et al. | ( |1999|) stream. 
b Subdwarf members of the newly disc overed kinematic groups. 
c Subdwarf members classified by |Klement et al. | (j2009|. 
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also discussed the possibility that such a stream is formed 
by the tidal debris of co Cen in the solar neighborhood, 
even though Navarrete et al. |(2015) have recently ruled 
out this hypothesis after detailed analysis of the chemical 
abundance of this group with respect to the well-known 
peculiar properties of u> Cen. 


4. SIMULATIONS 

We explore a simulated inner halo based on a set of four 
high-resolution numerical N-body simulations of minor 
mergers. We analyse the kinematics and orbital prop¬ 
erties of these simulations in order to investigate and 
characterise detectable signatures. 

It is useful to point out that these simulations are not 
an attempt to “fit” the observations, but they represent 
four merging events that we assume as representative, 
in terms of inclination and rotation, of the initial orbits 
of the satellites. In particular, these choices allow us to 
analyse the two cases suffering maximum and minimum 
dynamical friction. 


4.1. N-body simulations 

We use a set of high-resolution numerical N-body sim¬ 
ulations which simulate minor mergers of prograde and 
retrograde orbi ting satellite halos w ithin a dark mat¬ 


ter main halo (Murante et al. 2010). The main DM 


halo, which contains a stellar, rotati ng exponential disk, 


has a NFW radial density profile (Navarro, Frenk & 
White 1997), with a mass (M 2 qo = lO 1 ^ M 0 ), radius 


(A 200 = 105 kpc), and concentration (C 200 = 7.5), ap¬ 
propriate for a Milky Wav-like DM halo at redshift z = 0; 
the spin parameter is sefjjto A = 1. The satellite is rep¬ 
resented by a secondaryDM halo containin g a stella r 
bulge, and a Hernquist radial density profile (Hernquist 
1993); the spin parameter is set to A = 0. The mass ratio, 
AZprimary/-^satellite ~ 40, is similar to the estimated mass 
ratio of the Milky Way relative to the Large Magellanic 
Cloud. The main physical parameters of our simulated 
mergers are listed in Table [3] 

We consider prograde mergers, in which the satellite 
co-rotates with the spin of the disk, as well as retrograde 
mergers, with a counter-rotating satellite. We analyse 
two orbits: a low-inclination orbit with a 10° tilt with 
respect to the disk plane, and a high-inclination orbit 
with a 60° tilt. 

Initially, the particles of the small system (satel¬ 
lite galaxy) orbiting around the (otherwise static) disk 
galaxy are all strongly concentrated in space, and share 
essentially the same motion. The initial conditions (in¬ 
clination, position, and velocity) of the main system and 
the four impacting satellites, cosmologically motivated 
(Read et al. |2008[[Villalobos & Helmi ||2008), are sum¬ 
marised in Table |U 

From the grid of simulations by Read et al. (2008), 


we chose four impactors, all of which having the mass of 
the Large Magellanic Cloud. Larger masses would affect 
the stability of the stellar disk, and this is not consis¬ 
tent with a Milky Way-like galaxy. Conversely, smaller 
masses would produce minor signatures in our local halo 
sample. 


6 The cases A = 0 and A = 1 were both studied at lower reso¬ 
lution and the results compared; the differences are such that they 
have no bearing on the results presented in this paper. 


TABLE 3 

Physical properties of the Halos: Main system and orbiting 
Satellite. 


System 

M dm 

M, 

Ndm 

IV, 

ro 

T'd.isk 

Main 

10 12 

5.7 x 10 10 

10 6 

10 6 

4 

20 

Satellite 

2.4 x 10 10 

2.4 x 10 9 

1.1 x 10 5 

10 s 

0.709 



Note. — Column 1: Main galaxy/impacting system. Column 
2: DM mass, in Mg. Column 3: stellar mass for disk/bulge in the 
main/satellite, in Mg. Column 4: DM particles. Column 5: stellar 
particles for disk/bulge in main/satellite. Column 6: disk scale radius 
for the main halo, in kpc; Hernquist scale radius for the satellite, in 
kpc. Column 7: disk truncation radius, in kpc. 


The four simulations are compiled u sing t he public par¬ 
allel Treecode GADGET2 (Springe! 2005) on the clus¬ 
ter matrix at the CASPUR ( (Jonsorzio Interuniversitario 
per le Applicazioni del Supercalcolo) consortium, Rome. 
All systems were left to evolve for 4.63 Gyr (about 16 
dynamical timescales of the main halo). After this time, 
the four satellites have completed their merging with the 
primary halo. The final (x, z) distribution of the inner 
satellite star particles and host disk, in both the retro¬ 
grade and prograde cases, as well as for the high and low 
inclinations, is shown in Fig. [5j 


4.2. Dynamical friction and tidal stripping 

Any satellite can in principal be slowed by dynamical 
friction exerted on it by disk and halo particles. It is 
known that an object, such as a satellite, of mass M, 
moving through a homogeneous background of individ¬ 
ually much lighter particles wit h an isotropic veloc ity 
distribution suffers a drag force (Chandrasekar |[l943 ): 


F d = - 


4:TrG 2 M 2 pf(< v s ) In A 


where v s is the speed of the satellite with respect to the 
mean velocity of the field, Pf{< v s ) is the total density 
of the field particles wit h speeds less tha n v s , a nd In A is 
the Coulomb logarithm (Binney & Tremaine ||1987). 

We expect that, the higher the v s , the weaker is the 
dynamical friction force. Retrograde satellites are ex¬ 
pected to suffer weaker dynamical friction with respect 
to prograde ones, since in the first case the velocity of 
the satellite is opposite to that of the disk. As a conse¬ 
quence, prograde orbits decay faster. This effect is even 
more evident for low-inclination orbits. 

Another important effect that occurs during mergers is 
the tidal disruption of satellites. While tidal disruption is 
most important near the centre of the main halo, where 
the gravitational potential is changing more rapidly, dy¬ 
namical friction is exerted both by the main halo DM 
particles and by the disk star particles. 


4.3. Debris in the local halo 

We analyse the observational signature left by the 
satellite stars after selecting particles in a sphere of 3 kpc 
radius centered at the Sun (x = 8 kpc from the Galaxy 
center), and with |z| >1 kpc. This last constraint is 
introduced so that the simulated and observed samples 
can be compared within a similar volume of the inner 
halo. 
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Fig. 5. — Final configurations in the (x,z) plane of four minor merger events: depicted are the morphologies of the stellar distribution, 
i.e. the host disk (black) and the satellite bulge (colour), at the final time T = 4.63 Gyr of the simulations. Shown are the cases of 
low-inclination (10° tilt) retrograde (top left) / prograde (top right) orbit, and of high-inclination (60° tilt) retrograde (bottom left) / 
prograde (bottom right) orbit. The panels only display a randomly selected 10% subset of the total particles utilised. 


TABLE 4 

Initial conditions of the Main system and the four impacting Satellites. 


System Inclination 

Rotation 

(kpc) 

y 

(kpc) 

(kpc) 

V x 

(km s -1 ) 

Vy 

(km s 1 ) 

V z 

(km s -1 ) 

V 

(km s _1 ) 

Main 

0° 

— 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

0.00 

Satellite 1 

10° 

retrograde 

80.00 

0.27 

15.20 

6.30 

-62.50 

0.35 

62.82 

Satellite 2 

10° 

prograde 

80.00 

0.27 

15.20 

6.30 

62.50 

0.35 

62.82 

Satellite 3 

60° 

retrograde 

15.00 

0.12 

26.00 

1.20 

80.10 

2.00 

80.13 

Satellite 4 

60° 

prograde 

15.00 

0.12 

26.00 

1.20 

-80.10 

2.00 

80.13 


Note. — Inclination and rotation of the orbit, position and velocity components, and total velocity. 


Figure p] shows the kinematic distribution (velocity 
projections) of our simulated inner halo. The different 
colours indicate the association of the 3902 debris stars 
with different progenitors: the low/liigh-inclination 
retrograde satellites (761/616 green/red dots), and 
the high/low-inclination prograde satellites (966/1559 
blue/yellow dots). 

The angular momentum distribution of the satellite de¬ 
bris is shown in Fig. [7] Despite the chaotic build up of the 


parent halo, it appears that objects from accreted satel¬ 
lites remain confined in limited portions of the {L z , L xy ) 
plane. 

The satellite on a low-inclination prograde orbit (yel¬ 
low particles), which suffers more from dynamical fric¬ 
tion, quickly loses its orbital en ergy and proceeds to the 
inner regions of th e main halo (Byrd et al. ||1986 Mu-| 
rante et al. ||2010). For this reason less particles are left 


in the outer-halo/ see top-right panel of Fig. [5j 
It is worth noticing that the high inclination prograde 
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Fig. 6.— Kinematical (velocity space) distribution of the accreted component of the simulated Milky Way inner halo, i.e., 3902 particles 
in a spherical volume of radius 3 kpc centered on the “Sun” with \z\ > 1 kpc. Different colours indicate particles associated with different 
satellites: 60° retrograde/prograde (red/blue), 10° retrograde/prograde (green/yellow) colliding satellites. 
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Fig. 7.— Angular momentum distribution of the simulated Milky Way halo within 3 kpc of the “Sun”. As in Fig. [6] shown are the 3902 
particles accreted from four dwarf galaxies: 60° retrograde/prograde (red/blue), 10° retrograde/prograde (green/yellow) satellites after 
interaction with the simulated Milky Way. All of the marked regions have the same meaning as in Fig. [4] 
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satellite suffers the effect of dynamical friction as well, 
as a result of its co-rotation with the disk. This effect 
acts in producing a consistent mass of debris in the solar 
region having L xy ranging between 500 kpc km s _1 and 
1500 kpc km s -1 (blue points in Fig. [7|. 

On the other hand, retrograde satellites experience 
weaker dynamical friction and leave more particles in the 
outer-halo region, since their orbits have a longer decay 
ti me and longer periods. Thus, tidal stripping (see e.g., 


Colpi et al. 1999) can act longer and more efficiently 
when the satellite is still orbiting at high velocity, and 
we see that a better populated high-velocity tail results 
(compare Fig. [7 to Fig. 


10 ). 


The impact of dynamical friction on the two configu¬ 
rations considered for the retrograde satellites indicates 
that the high-inclination case is the one less affected by 
this force, which again results in efficient stripping when 
the satellite has high orbital velocity. Therefore, such 
stripping takes place over a large spatial region, and for 
the conservation of 6D phase-space density, by virtue of 
Liouville’s Theorem, we expect a small variance in veloc¬ 
ity space and in the plane of angular momenta. This is 
indeed observed for the red particles with respect to the 
green ones in Figs. [6] and [7] 

Finally, the effect of both gravitational feedback and 
dynamical friction on the satellites, which lead to loss 
of stars at different passages with different energies, is 
clearly evident for the case of low inclination prograde 
orbit in Fig. [ 7 ] at around L xy = 400 kpc km s -1 and 
L z = 1750 kpc km s -1 . 


5. “OBSERVED” SIMULATIONS 

Here we investigate the effects of observational errors 
on our simulated data, and show how more accurate kine¬ 
matic data to be provided by future surveys can improve 
detection and characterisation of halo streams. More¬ 
over, we compare actual observations with the distribu¬ 
tions of debris resulting from the four simulated satellites 
presented in the previous section, and discuss the orbital 
properties of the parent dwarf galaxies, possibly respon¬ 
sible for accreting on the Milky Way halo. 


5.1. Observational errors 

We perturb the original simulations by convolving the 
“true” data with two cases of error distributions. First, 
we adopt the accuracy of our SGKS catalogue as rep¬ 
resentative of the quality of current wide-held surveys. 
Then, we assume the mean accuracy expected from the 
forthcoming Gaia catalogue combined with complemen¬ 
tary deep spectrosco pic data from on-goin g and future 


surveys such as GES ( 

Gilmore et al. 

2012 

). 

The true positions and velocities o 

t eacJ 

a particle are 


first transformed into their astronomical observables (a, 
5, m — M , or 7r and /i a , fj,s , V r ); then the expected obser¬ 
vational errors are added to distance modulus (or directly 
to parallax, in the case of the Gaia-like simulation), ra¬ 
dial velocity, and proper motion, according to Table [5] 
The precision in distance is taken to a m -M = 0.4 mag 
(i.e., a d /d ~ 20%) for the photometric distances esti¬ 
mated from the SGKS catalogue, and to ov = 20 /ras 
for the final precision on trigonometric parallaxes mea¬ 
sured by Gaia. In proper motion, the precision is as¬ 
sumed to be 2 mas yr” 1 for ground-based observations, 


and 20 /ras yr -1 for Gaia. The precision in the radial 
velocity is taken to be 10 km s -1 for the SDSS measure¬ 
ments, and 1 km s -1 for the GES spectroscopic survey. 
These quantities are finally transformed back to observed 
positions vectors and space velocities. 

5.2. The inner halo model 

We explore a simulated inner halo based on the super¬ 
position of the four simulations of minor mergers and a 
smooth local component with the same kinematic prop¬ 
erties of the observed sample (Tabl e [Tj). 

C onsistently with t he fin dings of Helmi et al. (19991 


and Kepley et al. (2007), we assumed a debris total 


fraction of 10% within 4 kpc from the Sun. 

In the “true” (simulated) catalogue, of the 28 738 par¬ 
ticles with \z\ > 1 kpc, 24 836 are part of the local mock 
halo, while the remaining 3902 are debris from the satel¬ 
lites shown in Figs. [6] and [7] 

In the following discussion we focus on the particles 
of the accreted component that belong to the 10% high 
velocity tail. Table [6] reports the number of particles 
belonging to each satellite for the pure simulation and 
the other two cases accounting for observational errors. 

Figure [8]shows the region of the ({/, W) plane occupied 
by the 10% high velocity tail of the resulting simulated 
Milky Way inner halo (right panels), as superposition 
the accreted component (middle panels) and the smooth 
spheroid (left panels). The synthetic “observed” cata¬ 
logue shown in the top panels represent the current pic¬ 
ture, according to the SGKS error model. The bottom 
panels show the distribution of the high velocity particles 
as promised by Gaia. 

The upper panels indicate that distinguishing in ve¬ 
locity space the satellites that gave rise to each of the 
different moving groups with the extant data is a non¬ 
trivial task. On the other hand, as the inspection of the 
lower panels reveals, much of the substructures shown in 
the middle panels becomes visible again thanks to the 
superior precision that Gaia will achieve. 

5.3. Substructures in the correlation function 

In order to quantify the amount of kinematic substruc¬ 
tures present among the 2874 fastest-moving particles, 
we compute the cu mul ative velocity correlation function 
described in Sect. 3.2 The analysis is performed over 
three synthetic catalogues: the “true” simulation, and 
two lists derived from the true values after perturbing 
them with either SGKS-like errors or the errors expected 
for the Gaia/GES surveys. 

Figure [9] shows, using bins of width 5 km s -1 up to 
kinematical separations of 50 km s -1 , the results for the 
two-point correlation function £(v) for the “true” case 
(dots), the SGKS-like (squares) and Gaia/GES-like (di¬ 
amonds) catalogues. The clear signal in the first bins, 
peaking at Av ~ 15 km s -1 , evidences an excess of par¬ 
ticles moving with similar velocities with respect to what 
expected from a fully random sample. In the case of the 
pure simulation, the two-point velocity correlation func¬ 
tion attains a maximum signature of (£) = 0.68 ± 0.04. 
For the SGKS-like catalogue the maximum signal has an 
~ 80% drop to an “observed” value of (£) = 0.16 ± 0.04. 

The recovery of the intrinsic correlation signal is truly 
remarkable when looking at the correlation function for 
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Fig. 8.— The (t/, W) velocity distribution of the 10% high velocity tail of the simulated sample is shown in the right panels. The 
sample is limited to a spherical volume of radius 3 kpc located on the plane of the simulated Milky Way 8 kpc away from its center and 
with \z\ > 1 kpc. Of the 2874 particles in this sample, some are remnants of the four satellites accreted after 5 Gyr (middle panels), the 
remaining constitute the “background” of smooth inner halo stars (left panels). Different colours represent different progenitors as in Fig. pH 
Finally, the top panels were generated via convolution with current, i.e. SGKS, ground-based errors, while the bottom figures depict the 
results after convolution with the expected Gaia-like errors. 


TABLE 5 

Estimated/Expected errors for the SGKS and Gaia 

CATALOGUES. 


Catalogue 

distance 

proper motion 

radial velocity 

SGKS 

<*m-M = 0-4 mag 

2000 

10 

Gaia/GES 

cr 7 T = 20 fia,s 

20 

1 


Note. — Estimated errors (precision) in parallax (cr^, in 
/xas), distance modulus (cr m-M, in mag), proper motion (cr M , in 
fia,s yr —1 ), and radial velocity (cry r , in km s —1 ) for the SGKS and 
Gaia catalogues. 


TABLE 6 

The composition of the 10% High Velocity Tail of the simulated Milky Way Halo. 


Catalogue 

Halo+Debris 

Debris 

Satellite 1 

Satellite 2 

Satellite 3 

Satellite 4 

“True” Simulation 

2874 

1103 (0.38%) 

262 

201 

601 

39 

“Observed” SGKS 

2874 

835 (0.29%) 

233 

170 

406 

26 

“Observed” Gaia/GES 

2874 

1061 (0.37%) 

263 

191 

570 

37 
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Fig. 9.— The cumulative velocity correlation functions for the 
2874 halo particles shown in Fig. [8] The filled dots trace the corre¬ 
lation function of the pure simulation, while squares and diamonds 
depict the correlation after convolution with current (i.e., SGKS) 
or Gaia/GES-like errors, respectively. Error bars are derived from 
Poisson’s statistics of the counts. 


dots). Because of dynamical friction (cfr. Sect. 4.3), this 
satellite includes a low L xy component shown in the full 
sample_(Fig. [7|) that is not part of the high velocity tail 
(Fig. 10). Thus, these simu lated “ o bservations” suggest 
the possible presence in the Helmi et al. (1999) stream 
of debris with lower L xy yet to be discovered! 

Of particular interest is the case of the retrograde kine¬ 
matic groups. In fact, neither the high-inclination sim¬ 
ulated satellite nor the one at low-inclination appear to 
fairly match the observed Groups 3, 4, and 5, i.e. the 


black stars with L z < 0 in Fig. 10 
Actually, these three groups appear to occupy an inter¬ 
mediate region between the debris of the tw o sim ulated 
retrograde satellites. Furthermore, in Sect. |3.4 


we re¬ 


mark that the observed Grou ps 3, 4, an d 5 do no t we ll 


match t he stream s detected by Dinescu (2002) and Kep- 
ley et al. (120071). For this reason, we suggest that these 


three groups may represent the debris of an unique pro¬ 
genitor accreted along an initial retrograde orbit having 
an intermediate inclination in the range comprised be¬ 
tween 10° and 60°. In alternative, these groups could 
belong up to three different impacting satellites on ret¬ 
rograde orbits with inclinations in that same range. 

The results presented in this section show that the 
methodology proposed is certainly capable of detect¬ 
ing fossil signatures as kinematic substructures among 
high-velocity stars. From the data at our disposal, 
there is clear indication that more debris are found from 
dwarf galaxies on high-inclination prograde and retro¬ 
grade orbits, as well as on low-inclination retrograde 
ones. We have not identified any debris coming from low- 
inclination prograde satellites and this might be a limita¬ 
tion intrinsic to the methodology of looking at structures 
in the space motions of very high velocity stars. Future 
work will have to investigate these issues. 


the Gaia-like case: in fact, we measure (£) = 0.46 ±0.04, 
corresponding to 68% of the original signal. Figure [8] 
provides a nice “visual” confirmation of the recovery in 
substructure visibility. 

6. ON THE NATURE OF THE HIGH VELOCITY 
DEBRIS 

As the space of adiabatic invariants is important 
to gain more insight into the pr oper ties of the kine¬ 
matic substructures detected (Sect. 3.4), we compare the 
(. L z ,L xy ) distributions of the observed groups with the 
results of the simulations, taking into account the effect 
of the observational errors. This is shown in Fig. |10[ 
where the top panel corresponds to the SGKS-error sim¬ 
ulation, while the bottom panel reproduces what will 
hopefully be seen with the final Gaia catalogue. 

The black star symbols in the upper panel of Fig. [10] 
represent the 67 high velocity objects we found from 
our statistical analysis in the same volume and shown 
in Fig. [4] as colored stars. With current data, different 
satellites mix over some regions so that a discrete 
classification is not always straightforward. The bottom 
panel of Fig. [10] clearly shows that this situation is 
highly improved with Gaia-like data. 

We see t hat o ur Groups 1 and 2, corresponding to the 


7. CONCLUSIONS 

We have explored the Solar neighborhood of the Milky 
Way through the use of spectro-photometric data from 
the Sloan Digital Sky Survey and high-quality proper mo¬ 
tions derived from multi-epoch positions extracted from 
the Guide Star Catalogue II database. A sample with 
accurate distances, space velocities, and metallicities is 
selected as a tracer of the inner-halo population resulting 
in 2417 subdwarfs with [Fe/H] < —1.5 and |z| > 1 kpc 
within 3 kpc of the Sun. This set is then analysed to iden¬ 
tify and characterise kinematic streams possibly arising 
from merging events. 

We have found statistical evidence of substructures in 
the space motions of the 10% fastest stars, confirming 
the existence of 5 moving groups. 

In angular momenta space, the two prograde groups 
we have identified (Groups 1 and 2 in Table pi) appear 
confined i n the region enc o mpass ing the stream nrst iden¬ 


tified by Helmi et al. (1999) among red giants and 


RR Lyrae within 1 kpc of the sun. Our analysis found 
25 additional subdwarf members belonging to that_same 
stream: 4 are in common with those found by Klement 
et al. (2009), while the other 21 are newly discovered 
members. 

Of the remaining groups, the most counter-rotating 
one (Group 4) partially ove rlaps with the region of “ret- 


stream of Helmi et al. (1999), are consistently asso- rograde outliers” found by Kepley et al. (2007), while 


dated with the high inclination prograde satellite (blue a dozen stars belonging to Group 4 and Group 5 fall in 
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-4000 -2000 0 2000 4000 

L z (kpc km s" 1 ) 

Fig. 10.— The 10% high velocity tail component belonging solely to the 4 satellites; convolved with current ground-based errors (top, 
835 particles) and with the expected Gaia errors (bottom, 1061 particles). Black star symbols are the 67 fast moving debris stars uncovered 
from the analysis of the SGKS sample. Solid and dashed contours have the same meaning as in Figs. [j]and[7]. 
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the regio n of the mildly retrograde stream detected by 
Dinescu (2002]). 


Comparison to our high resolution N-body simulations 
confirms that the two groups associated with the Helmi 
stream are likely fossil remnants of a dwarf galaxy which 
co-rotates with the disk of the Galaxy and moves on a 
high-inclination orbit. 

As for the three retrograde groups (3, 4 and 5 in Ta¬ 
ble [2]) , they may be debris of an unique progenitor ac¬ 
creted along an initial retrograde orbit having an inter¬ 
mediate inclination in the range 10° -j- 60°. However, we 
cannot exclude that these groups belong to different im¬ 
pacting satellites on retrograde orbits with inclinations 
within that same range. A more detailed analysis of the 
chemical abundances of the three detected groups, as well 
as more quantitative comparisons to extended simula¬ 
tions, are necessary to resolve this issue. 

In any event, the fastest objects appear with pos¬ 
itive and negative L z values (i.e. prograde and ret¬ 
rograde motions, respectively) in the angular momen¬ 
tum ( L z ,L xy ) regions for high-inclination orbits ( L xy > 
1500 kpc km s -1 ). On the other hand, for low- 
inclination both observations and simulations show that 
the fastest objects appear only on retrograde orbits (e.g., 
Fig. 10 top panel). This asymmetric distribution is sug¬ 
gestive of the role played by dynamical friction during 
accretion. 

In anticipation of the much improved data expected 
over the coming years, in particular the Gaia catalogue 
and the new ground-based spectroscopic surveys, we also 
investigated the impact of observational errors in our dy¬ 
namical simulations. The analysis indicates that (see the 


relevant panels of Figs. [8] and 10) Gaia will greatly influ¬ 
ence these studies: for, velocity and angular momentum 
distribution will be almost completely dominated by the 
physics we are trying to recover, i.e., the dynamical his¬ 
tory of the merging events. 

At that point, full grids (in, e.g., inclination and 
amount of rotation) of prograde and retrograde high res¬ 
olution satellite simulations will be required to precisely 
characterise the debris detected. Then, we will be able 
to number the merging events for direct comparison with 
the predictions of the (A)CDM theory and its associated 
merging paradigm. 

In conclusion, the results shown might lead us to 
claim that the inner halo might have “seen” only two 
events; however the large uncertainties in the extant 
data, mostly observational, do not exclude the possibility 
that the events might be as many as four, and perhaps 
more given the intrinsic difficulty of our technique to deal 
with low inclination prograde mergers. 
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Academy of Sciences through 2015 CAS President’s In¬ 
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